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The  following  two  multiple  integrals  are  evaluated 
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Also  the  asymptotic  behavior  of  the  two  integrals  as  |x 


and  as  jx|— «>  is  given. 


MULTIPLE  INTEGRALS  INVOLVING  PRODUCT 
OF  MODIFIED  BESSEL  FUNCTIONS  OF  THE 
SECOND  KIND 
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1.  Introductory.  It  is  proposed  to  establish  the  two  following  integrals: 


If  m  =  2,  3,  4, 
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where  x  is  real  and  positive  and  the  symbol  E  means  that  in  the  expression 

following  it  i  is  to  be  replaced  by  -i  and  the  two  expressions  are  to  be  added. 
If  R(  k_  =  ±  jx)  >0,  r  =  1,  2,  . . . ,  m-1  , 
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where  m  =  2,  3,  4,  . . . ,  x  is  real  and  positive  and  the  symbol  7)  has  the 

same  meaning  as  before.  The  function  appearing  in  { 1)  and  (  2)  is  Mac  Robert’ s 
E-function,  the  definitions  and  properties  of  which  are  to  be  found  in  [  1]  ; 
pp.  348-352  and  in  [2]  ;  pp.  203-206.  The  following  formulae  will  be  required 
in  the  proofs: 

If  R(s±p.)  >0,  then  {[3],  p.  197  and  p.  ?) 


/  xs  1  k  ( x)  dx  =  2s  2  r^)  n^) , 
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If  P  >  q  +  1,  then([l],  p.  409) 
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where  jarg  z  j  <~  (p -q-J-1)  tt;  the  prime  in  the  product  sign  signifies  that  the 

Ml 

factor  for  which  s  =  r  is  left  and  the  asterisk  in  the  F  function  means  that  the 
parameter  or  -  +  1  is  omitted. 

If  P  <  q,  then  ([1],  p.  352) 
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If  m  is  a  positive  integer 
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1  i  1 

r(z)  r(z+“)  ...  r(z+— )  =  (2tt)  m  r(mz) 
m  m 


(V 


If  m  is  a  positive  integer,  then  ([4],  p.  759) 
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where  p  >  q  4-  1  and  | arg  z |  <~( p-q-1)  ir  . 

If  iarg  z|  <  ^(  p-q4-l)  v,  then,  ([1],  p.  374) 
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where  the  integral  is  taken  up  the  pL-axis,  with  loops,  if  necessary,  to  ensure 
that  the  pole  at  the  origin  lies  on  the  left  and  the  poles  at  a^,  . . . ,  a  to 

the  right  of  the  contour.  The  contour  must  be  modified  if  p  <  v  4  1  and  p  -  q  +  1 
|z|  <1  . 

Also  the  proof  depends  on  the  expression  in  terms  of  E-function  of  the  gen¬ 


eralized  E-function 
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where  i  and  m  are  positive  integers  and  the  contour  passes  up  the  Tj-axis 
from  -oo  to  +co?  with  loops,  if  necessary,  to  ensure  that  the  poles  of  the  in¬ 
tegrand  at  the  origin  and  at  p  .  -1,  p  -1,  ....  p  -l  lie  to  the  left  and  the 
poles  at  a^}  a . . . ,  to  the  right  of  the  contour.  When  necessary  the  con¬ 
tour  is  bent  to  the  left  or  to  the  right  at  both  ends  until  it  is  parallel  to  the 
£-axis. 


This  expansion  {[1],  p.  419)  is 
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where  the  prime  and  the  asterisk  denote  that  the  factor  sin(p  -p  and 

q-*-s  q+s 

the  parameter  p^+s~p^+s+l  are  omitted  and  ca  is  equal  to  1  or  e  according 
as  I  +  m  is  even  or  odd. 

To  familiarize  ourselves  with  the  E-function  the  following  relations  may  be 
worth  noting;  they  will,  however,  not  be  used  here.  From  the  definitions  ( 5) 
and  ( 6)  of  the  E-function  it  is  clear  that  the  E-function  is  immediately  related 
to  the  generalized  hypergeometric  function  and  reduces  to  simple  expressions  in 
the  ordinary  or  Gauss  hypergeometric  function  when  p  =  2,  q  =  1  .  For  p  =  1 
and  q  =  1  it  is  also  evident  that  the  E-function  reduces  to  the  confluent  hyper¬ 


geometric  function  or  Kummer 1  s  function.  The  case  p  =  2,  q  =  0  yields  the 
relations  [see  [1],  p.  351]. 

cosmr£(^  +  n,  ^-n!!  2z)  =  *s/(  2ttz)  eZK  (  z) ,  (12) 

z  z  n  j 

—  2 

E(“  -  k  +  m,  “r-k-m  !!  z)  =r(“-k-m)  r(Tr-k  +  m)  z  W  ( z)  (13) 

z  z  z  z  k,  m 


where  K  ( z)  and  W.  (  z)  are  the  modified  Bessel  function  of  the  second  kind 
v  k,m 

and  Whittaker  function.  Also  it  is  immediate  from  the  definition  of  the  E-function 

-l/z 

that  for  p  =  q  =  0  the  E~function  is  just  e  and  we  have 

_1_ 

E( !  I  z)  =  e  Z  (14) 


More  complicated  parameters  in  the  E-functions  lead  to  the  equivalence  of  the 
E-function  with  products  of  Hankel  functions,  with  Lommei  functions  and  with 
products  of  Whittaker  functions.  Some  examples  of  this  are: 
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§2.  Proofs  of  the  integrals. 

To  prove  ( 1) ,  apply  (  4)  to  replace  K 
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1  m-1 


on  the  left  of  ( 1)  by 
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Then,  on  changing  the  order  of  integration  and  evaluating  integrals  by  means  of 
(3) ,  the  multiple  integral  becomes 
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where  we  have  used  the  relation 


m)  rd-e)  = 
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sin 


Now  apply  (  9)  with  p  =  2m  +  1,  q  =  0  and  so  obtain  ( 1)  . 

Again  on  substituting  for  K  (  x\  . . .  X  )  on  the  left  of  (  2)  from  ( 4)  and 
changing  the  order  of  integration  as  before,  the  multipie  integral  is  found  to  be 
equal  to 
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and  on  applying  (11)  with  m  =  2,  q  =i  =  0,  the  expression  on  the  right  of  (  2) 


is  obtained. 


§3.  Particular  cases. 


In  (1)  take  k.  =  -^-,  u.  =  p  (r  =  1,  2, . . m-1)  and  it  becomes 
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Formula  (  21)  is  a  generalization  of  Hardy' s  formula  { m  =  2)  ( [  5] ,  p. 


namely 


/k  (M  K  (*)dk  =  *K,  (2r(x))  , 


(  22) 


where  x  >  0  . 
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In  { 1)  take  k  =  —  ,  n  =  v  =  ±  \  ( r  =  1,  2,  . . . ,  m-1) ,  apply  the  formula 
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and  it  becomes 


m-1  oo  —  - 1 

„  r  r  .  m 
n  J  e  \ 

r=l  0 


dX  e 
r 


X_  *  •  •  X  . 
1  m-1 


=  2 


~2  _m 


1  1 

2 m  2 
x 


v  J_  y  iFli  i  1,1  1 ,  i 

2ir  i  4  3  43  4  2m  ’  *  2  ’ 


=  x 


1 

2  1 


2t>  H  . 


1  ,  m-1  1  ,  m-1  ...  iTT  x  . 

.....  .  ..  j  *  *  *  e  .  ) 

3  4  m  4  m  '  2m 

V  1  w  1  1.1  1  ,  m-1  ...  i"  . 

H  .  i  E("2  j  ~2  +  m>  *  *'J  "2  *  m  }  1  *  *  6  X) 


i,  -i 


=  xr  E  ~  E(-l,  -It—  ,  ..  .,  -1  -r— ,  1  ::  e17r  x) 
2it  y  .  i  m  3  m 


i,  "I 


1 


1  ,  1 
—  (m-1)  m 

,  „  ,  2  2  ...  u  m  , 

=  x  (  2 it)  m  E(  -'m,  1  . .  e  x  m) 


by  ( 8)  with  <*  =  -1,  £*■  =  — ,  q  =  0  ,  and  so  we  obtain 
J  1  ’  2  m  ’ 


.  1  i  „  m-1  , 

oo  —  - 1  oo  —  x.  — -  -  1 

r  1  ,m  „  r  m-i  .  m  .. 

/  e  X,  dX,  —  /  e  X  dX  .  e 

1  1  m-1  m-1 

1  1  — 

,,  .2(m'1)  ~Z  -mxm 
=  (  2w)  m  e 


i 

i  m-l 


(24) 


where  x  is  real  and  positive. 

Again  in  ( 1) ,  take  m  =  2,  and  expand  each  E- function  on  the  right  of  ( 1) 
by  means  of  (5)  and  combine  the  two  resulting  expressions,  so  getting, 
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where  x  >  cr  and  each  of  the  symbols  Z  ,  Z  has  the 

|i,  -|i  V,-V 


same 


meaning  as  before. 


In  (  2)  j  take  m  =  2,  and  get 


CO 

/  Xk_1  K,(  X)  K  (xX)dX  =  2k~3 


Z-p-  k+v-n .  .  k-v-LL 

x  n—r^)  rfn) 


-H- 


F  k-v-p.  2. 

2  1  2  3  2  J  1  "  *  ;  x  ) 


Now  sum  the  two  hypergeometric  series  of  Gauss  by  the  formula 
F(  a,  b;c; 2)  -  Zu-z, -  ^(a,  c-b,  a-b  +  1;  £  , 

where  |arg(  1-z)  |  <  tt, 

and  so  obtain  Titchmarch1  s  formula  ( [6]  p.  98)  namely 


(  25) 
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f  xk  1  K  (X)  K  (xX)dX  =  2k"3  {r(k)  }-1  x  V~k 
0  v  ll 


r{M|±V}  r{Mp}  r(k=p} 

p  fk+E±Z  Mzzii.  v  .  !  _  _L  ^ 
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x 

where  R(  k  ±  p.  ±  v)  >  0  and  x  is  real  and  positive. 

In  (  2)  take  m  =  3,  and  obtain  the  double  integral 
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where 

R(  k.  ±  v  ±  p.)  >  0,  R(  k  ±  v  ±  p)  >  0?  x  is  real  and  positive. 

1  1  Lt  L> 


(26) 


(27) 


§3.  Asymptotic  behavior  of  the  integrals: 

Now  we  are  in  a  position  to  find  the  asymptotic  expansion  of  the  two  multiple 
integrals  (1)  and  (  2)  for  large  and  small  values  of  |x| ;  for  the  asymptotic 
expansion  of  the  E-function  is  given  by  Mac  Robert  ( [1] ,  p.  358)  . 

If  p  >  q  +  1 


E(p;  a  :  q;  pg  :  z) 


the  first  m  terms  of  the  divergent  series 


r( aj  ...  r(ap)  l 
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(28) 
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where  R  ,  is  a  remainder  which  tends  to  zero  as  j  z  |  —  «  and  m  —  °o  . 

m-rl 

And  the  other  functions  that  occur  can  be  expressed  in  terms  of  ordinary 
generalized  hypergeometric  functions  whose  asymptotic  expansions  have  been 
investigated  by  several  writers  (Barnes,  [7]  and  [8],  Wright  [9]).  However, 
an  alternative  proof  for  the  asymptotic  expansion  of  the  integral  (1)  as  |x|— «> 
will  be  established.  Thus  we  shall  prove,  for  large  values  of  |x|,  the  follow¬ 
ing  asymptotic  expansions 
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where  the  coefficients  and  M2  do  not  depend  on  x  but  are  complicated 
functions  of  the  k^,  v  p  ( r  =  1,  2,  . . . ,  m-1)  , 
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where,  for  simplicity,  we  have  taken  k  =  k(  r  =  1,  2,  . . . ,  m  -  1,  (  30} 

To  prove  {  29) ,  substitute  from  ( 9)  for  the  E-functions  appearing  on  the  right, 
then  it  becomes 
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We  next  show  that 
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To  see  this  we  substitute  the  value  of  S{  s)  in  ( 32) ,  the  right  hand  side 
becomes 
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which  is,  in  virtue  of  (31) ,  equal  to 
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We  now  use  (7)  with  2m  instead  of  m  to  obtain 
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With  (32)  established,  we  utilize  the  results  of  Barnes  on  the  asymptotic  ex¬ 
pansions  of  the  generalized  hypergeometric  functions  (see  [7],  pp.  296-297, 
and  [8],  pp.  80,  108,  110)-  for  the  asymptotic  expansion  for  the  expression  on 
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the  right  of  (  32) ,  which  holds  when  x  is  large  and  so 
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and  so  (29)  is  established. 

By  applying  (  5) ,  formula  (30)  can  be  established  for  large  values 
Of  course  the  leading  term  in  the  asymptotic  expansion  of  (  29)  is 
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We  now  consider  the  special  case  k  =  —  and  v  =  u.(  r  =  1.  2, . . 
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and  show  that  the  result  (34)  agrees  with  the  asymptotic  expansion  of 
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which  agrees  with  the  asumptotic  expansion  of  the  right  hand  side  of  (  21) 
because 
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Again  we  can  write  the  asymptotic  expansions  of  the  integrals  ( 1)  and  (  2) 
as  1  x |  - *■  °°  .  Tor  the  asymptotic  expansion  of  ( 1)  as  |x|  — co  is  obtained  by 
expanding  each  E-function  on  the  right  by  means  of  (5) ,  and  combining  the 
two  resulting  expressions  by  factoring  out  common  terms. 

The  asymptotic  expansion  of  the  integral  (  2)  can  be  written  at  once  in 
virtue  of  (28)  in  the  form 
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